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Abstract 

We present an analysis of radiation induced instabilities in an expanding plasma with consid- 
erable presence of dust particles and equilibrium cooling. We have shown that the equilibrium 
expansion and cooling destabilize the radiation condensation modes and the presence of dust par- 
ticles enhances this effect. We have examined our results in the context of ionized, dusty-plasma 
environments such as those found in planetary nebulae (PNe). We show that due to the non-static 
equilibrium and finite equilibrium cooling, small-scale localized structures formed out of thermal 
instability, become transient, which agrees with the observational results. The dust-charge fluctua- 
tion is found to heavily suppress these instabilities, though in view of non- availability of convincing 
experimental data, a definitive conclusion could not be made. 
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I. INTRODUCTION 



Thermal instability is long thought to be a relevant astrophysical process which is respon- 
sible for existence of smaller, non-gravitational condensations such as those found in solar 
prominences, interstellar clouds, and planetary nebulae (PNe). Earlier works on thermal 
instability include those of Parker^, Weymann^, and Field^, among which. Field's seminal 
work on thermal instability is considered to be a comprehensive treatise on the subject till 
date. Recently, Gomez-Pelaez and Moreno-Insertis^ have addressed the process of thermal 
instability in a gravitating medium undergoing expansion. Several other authors have con- 
sidered occurrence of thermal instability in different contexts and parameter regimes. An 
analysis of the thermal instability of an optically thin plasma in the nonlinear regime is 
considered by Steele et. al.- Some of the very recent works involving numerical simulations, 
address the issues of dense structure formation in a hot protogalactic environment^"^ and the 
interplay between thermal and magnetorotational instability (MRI) in interstellar mediai^>^. 

It is also now well known that the dust particles constitute an ubiquitous and important 
component of many astrophysical plasmas including interstellar clouds, stellar and planetary 
atmospheres, planetary nebulae, giant H-II regions (GHR), hot interstellar matter (ISM), 
etc. The presence of dust in an ionized astrophysical structure, besides being responsible for 
introduction of new wave modes related to the dust dynamics, can also significantly modify 
the thermal structure of the plasmap^iii. One of the important ways dust can modify the 
stability properties is by capturing energetic electrons and cooling through radiation^i^ in 
the infrared region and we expect that the presence of dust particles modify the correspond- 
ing growth rates significantly. The charge fluctuation dynamics of the dust particles is also 
found to introduce new acoustic modes in the plasmeJ^^. Besides, charging of the dust 
grains makes them an integral part of the multi-component plasma which also helps main- 
tain the colloidal nature of the plasma. Apart from the physical influence of the presence of 
dust particles in such astrophysical environments, the scenario becomes more complicated 
when one considers the dynamic behavior of the background. For example, the planetary 
nebulae are known to have severe equilibrium expansions, in many cases with supersonic 
expansion velocities^'i^,. It is therefore of particular interest to review the process of con- 
densations in ionized plasmas taking into consideration the effect of the dust dynamics and 
their participation in the overall stability with a non-static background. In this work, we 
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have undertaken a systematic study of the thermal condensation process in such an ionized 
plasma environment which undergoes equilibrium expansion and cooling having a consider- 
able presence of dust particles which constitute another charged species of the plasma. We 
would like to view the outcome of this analyses in the context of condensations observed in 
planetary nebulae (PNe). 

The problem of existence of dust particles inside the hot ionized environment of a PNe 
is yet to be understood fuUy^^. It has however been proposed that majority of the dust 
grains could be destroyed in the hot ionized regions of plasma in a PNe, as it evolves, either 
because of spallation by hard UV photons or shocks. Subsequently, a PNe or a GHR should 
be free of dust grainsi^i^S. However, in many cases the destruction of dust particles inside 
an ionized environment is not complete. For example, recent observations of the infrared 
spectra of the evolved nebula NGC 6445^1^, it has been shown that only a little destruction 
of the dust grains inside the ionized region of the nebula could have occurred and dust-grain 
separation inside the ionized region is not plausible^. Even earlier attempts to detect Ca 
lines in PNes suggest that Ca is probably locked in the dust grains present inside the ionized 
environmental^. Theoretical investigations, so far, indicate that dust particles are gradually 
accelerated to higher and higher velocities by the radiation pressure and are being expelled 
from a nebula^^i^i^ which however do not take into account the mechanism of charging of 
the dust grains. The charging of the dust grains freezes them into the ionized gas^^, about 
which we have already referred. 

In Section 2, we describe our model equations with background cooling and expansion. 
In Section 3, we apply the linear perturbation theory to examine the thermal condensation 
modes. In Section 4, we consider the case of static dust charge with the help of a WKB 
formalism and derive the linear dispersion relation. We numerically examine the dust-charge 
fluctuation dynamics in Section 5. Finally, we summarise our observations in Section 6. 

II. MODEL EQUATIONS 

Below, we write down the equations for a weakly coUisional, fully ionized dusty plasma 
with an expanding background. The equations are namely, the equations of continuity, 
momentum, and energy conservation. 




(1) 
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rrijUj + '^j ■ ^'^i^ = ^'^Pj ~ e%V0 — mjUj'Vtlj, (2) 

3 dT ' 

where the subscript j = e and d for the electrons and dust particles, Xj is the thermal 
conductivity of the species, and and are the electrostatic and gravitational potentials 
defined by their respective Poisson's equations, 

= -47r[e(n, - n^) + Qna], (4) 
VV = ^irGmand, (5) 

Q being the total charge on the dust particles. The ratio of specific heats 7 = 5/3 in these 
equations and other symbols have their usual meanings. As mentioned before, we have 
assumed the massive dust particle to be negatively charged. The second term on the right 
hand side of Eq.Q represents the net heat-loss through radiation^. Note that, in writing 
these equations, we have assumed that the ions are Boltzmanian, 

Hi = noe~~- . (6) 

This assumption of thermalization of ions comes from the fact that at temperature above 
~ 10^ °K, it is effectively the electron energy, which gets radiated and the radiative cooling 
of the electrons prevents them from attaining thermal equilibriumi^i^^. Besides, as we are 



interested in the region of perturbation with phase velocity of the order of ~ JTe/md] the ion 



thermal velocity (~ \jTe/mi) can be assumed to be much higher, so that any temperature 
fluctuation in the ions are assumed to be quickly equilibrated. The pressures for the electrons 
and dust particles are given by pj = njTj. In what follows, we however neglect self-gravity, 
as the regimes where our analysis might be applicable are not sufficiently massive to be 
affected by self-gravity—. 



A. Background (equilibrium) expansion with net cooling 

We consider now a homogeneous and uniformly expanding background and assume that 
self-gravity is too weak to affect the expanding plasma so that the equilibrium is character- 
ized by an expansion parameter e{t}^^, 

r = e{t)x, e{t) = const. (7) 
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where x is the Lagrangian coordinate. Due to this non-static background, the equihbrium 
quantities naturally become time dependent. Note that when self-gravity is not negligible, 
the rate of expansion is not constant and is given by, 

4 

where po is equilibrium matter density. The equilibrium background velocity Vq is given by 

Vo{r,t) = -j. (9) 

This expanding equilibrium is mainly characterized by the equilibrium density rtjQ and tem- 
perature Tjo, 



dt e 



3-njo, (10) 



|.ogr,„^-2(i4„,g. (u) 

where the equilibrium quantities are denoted by a subscript '0'. In the expression for evo- 
lution of equilibrium temperature Eq. (fTT|) . there are essentially two effects causing temper- 
ature to change — equilibrium background expansion [the first term in Eq.l fTTll ] and net 
equilibrium cooling (the second term). 



III. THE LINEAR PERTURBATION ANALYSIS 



We perturb the equilibrium with a small electrostatic perturbation and for any arbitrary 
physical quantity f{x,t), we write the corresponding perturbation as 

A(a;,t)~ A(t)e'^-^ (12) 

where k is the Lagrangian wave vector. In doing so, we have retained the time dependent part 
of the perturbation exclusively with the amplitude and Fourier decomposed the space part^. 
We spht the perturbed velocity into two components, namely parallel and perpendicular to 

k 

Vji{t) = -^v.ji{t) + Vji^{t), (13) 

where Vjiit) = k ■ Vji. 
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From the electron and dust continuity equations, the perturbed densities can be written 



as, 



dfij 

— = -^kv,„ 



(14) 

where the quantities with a ' " ' are perturbed quantities normahzed by their equihbrium 
values and k = k/e. If we neglect the inertia term of the electrons in the momentum 
equation, Eq.([2]), compared to the massive dust particles, the perturbed electron pressure is 
just the perturbed electrostatic potential. 



Pel 



e(pi 

TeO 



(15) 



in which, we have normalized the potential with the equilibrium electron temperature. Note 
that the temperature is expressed in energy units. As the ions are assumed to be Boltzmanian 
[Eq.([n])], the perturbed ion density can simply be expressed as. 



(16) 



where r = Teo/Tio is the ratio of the electron and ion equilibrium temperatures. 
From the electron energy equation Eq.([3]), we have. 



3 dTe „ 



PeiXok^TeO + meUeoC-eO + meUeoTj^Cn 

-heixok^Teo + meUeoT^oCp^ 



0, 



(17) 



where the ()i denotes the first order perturbation of the respective quantity and the equi- 
librium derivatives of the radiative loss function are defined as. 



r 



Pe 



d 




dTeo 


\t7o 


d 




dTeo 


\t7o 



The quantities inside the ()i are given as, 

(PeV -V^^ = Peo(V -Vel) - p, 
dTei 



3 dTeo „ 

■KneO—— + meUeoL-eO 



Hp 



dTe 
dt 



neo- 



+ neo{Pe-Te 



dt 
, dTeo 



(18) 
(19) 

(20) 
(21) 



dt '^"^"^ ^' dt 

Finally using the perturbed electron continuity equation, Eq (|T^ . we can write Eq. (|T71) for 
the perturbed electron temperature as. 



dpe 
dt 



.5- 



i::kvei - Pe{r ^ + uJn) + rieir ^ + Up), 



(22) 



where r-^, ci;„^, and ^ are the characteristic time-scales for the thermal conduction by the 
electrons, isochoric, and isobaric electron perturbations^, given by 



^n,p — 7T^eT'e0^jie,pe- (23) 



^ 3 ' 3- 
The quasi-neutrality condition with negatively charged dust particles is given by, 

he6e -ni + {hd + Zd)Sd = ^, (24) 

where we have taken the possibility of finite dust-charge fluctuation. The quantities 5e = 
neo/nio and 6d = ndoZdo/riio are the ratios of equilibrium electron and dust densities to the 
ion densities and Q = qZd is the dust charge with Zd being the dust charge number. The 
dust charging equation is given by^, 

^=/ei + /a, (25) 

where, Ijis are perturbed electron and ion currents. The above charging equation can be 
ultimately reduced to, 

—1— + vQdi = \Ieo\ , 26) 

at Kriio rieoJ 

where rj is natural decay rate of charge fiuctuationsi^. 

Finally, we write the dust momentum equation Eq.([2]) as, 

= --Vdi + iZdoUJdCd(t)i, (27) 
at e 

where 

Ud = kcd (28) 



is the frequency of acoustic perturbation, and q = yTeo/rud. By virtue of the non-static 
equilibrium, all the ujs i.e. ujn^p^ have become time dependent quantities. 

The complete closed set of equations for linear perturbation of our model is given by 
Eqs.([Tl-(Il6]), (122D, ([26D, and ^ 



A. Analysis of the perturbed system 

As the linearised differential equations become non-autonomous due to equilibrium ex- 
pansion and net cooling, the problem, in general, is not amenable to standard Fourier decom- 
position for obtaining the dispersion relation, which requires rigorous numerical treatment as 
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exact analytical solutions do not exist. However, depending on the time-scales of the prob- 
lem, we can apply approximation techniques viz. WKB approximation:^»^ to gain valuable 
insight into the nature of the problem. Recently Nejad-Asghar and Ghanbari^ have con- 
sidered the problem of non-static equilibrium for radiation instability in molecular clouds, 
where they have used an ansatz of exponential form for the normalized perturbed quantities. 
In our notation, their formalism can be expressed as. 



where the quantity uj is analogous to normal mode frequency, which is decoupled from t. 
However, as the equilibrium is time dependent, the normal modes of the system become time 
dependent too, so that uj = uit) and a decomposition of the kind expressed in Eq. (!29!) is not 
possible. As we show in the appendix, perturbation amplitude of densities based on Eq.( l29|l 
does not lead too far and grossly disagrees with the numerical analysis and the validity of 
the corresponding results are very limited. 

In the following Section, we use the WKB approximation to analyse the perturbed system 
in the limit of static dust-charge. 

IV. STATIC DUST CHARGE : WKB APPROXIMATION 

In order to simplify the mathematical complexity of the problem, we first consider the 
case of static dust charge. In this limit, we can reduce the relevant equations of Sec. II to 
form a third order differential equation in any of the perturbed variables, which can be 
conveniently expressed as. 



rid ~ 




(29) 



+ ai(t) 



dfid 



+ aQ{t)hd = 0, 



(30) 



dt^ 



dt 



where the coefficients aj(t)s are given by. 



a3(t) = 35e + 5r(t) 



(31) 
(32) 
(33) 
(34) 



a^it) = 3(5ef^n + 3r(t)f]p + 65et^c + 2cjj95e + 10r(t)], 
ai(t) = 2a2(t) - LJ,[96, + lOr(t)] + 56dZdo{t)ujl 



ao{t) = 36dZdo{t)u, 



d^'^pi 
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Various frequencies used in Eq. (l30|) are defined as 



Qn = ujn{t) + r-\t), (35) 
fip = ujp{t) + T-\t). (36) 



The expansion and cooling time-scales are given by u^l, 



u^s = [-) , = (^-m.^ j . (37) 

Different time-scales of the system are dictated by ci;^^p^c,e,x' ^ clear separation of which is 
essential for WKB approximation. 



A. Linear dispersion relation 

Going by the standard WKB theory^^, we expand the eigenfunctions nd{t) of Eq. flSU]) as, 



hd(t) ~ exp 



1 oo 

TE'^"^n(t) 



5 



6^0, (38) 



where 5 is a small parameter and Si{t) are arbitrary functions of time. We break up the 
coefficients ai{t)s in time as, 

a,(t) = + a^^e^' + af ^e^^ _^ . . . ^ (3g) 

where e is the smallest time-scale of the system. We assume a-priori that the different 
coefficients aj(t)s can be ordered as, 

as{t) = , (40) 

a,{t) = 4°)+4^)e-\ (41) 

a,{t) = af^+ai'^e-i + af^e-^^ (42) 

ao(t) = + a[,'^e-i + aj'^e"' + a^^^e-\ (43) 

Assuming e ~ 5 = e, we substitute the above approximations along with expansion given by 
Eq. (1391) in Eq. flSU]) . After equating the like powers of e in the resultant relation, we obtain 
the linear dispersion relation for the system, from Eq. fl30l) as. 



"3 



N-1 N N-j 

3 Y: S',S';,_,_^ + S%_, + S'^S[S'^_i_^ | = O, iV = O, l, . . . , oc<44) 

j=0 i=0 1=0 



'N-l \ I ( ^ 

J=o J \j=o 
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where 



af\ Sj = 0, for j < 



(45) 



i = 0,1,2,3 

We note that apart from the basic assumption of well separated time-scales of our problem, 
two essential conditions for validity of the WKB expansion are, as e ^ 0, 



We express the WKB function Sn(t) as an integral of an unknown function ipnit), 



where tq is the shortest time-scale associated with the physical system. The WKB dispersion 
relation, Eq. fH^ can now be examined in different regimes depending on their characteristic 
time-scales. 

B. The dust-acoustic domain (wc ~ ~ „ ~ <C uJd) 

In the dust-acoustic domain, we assume that the characteristic time-scale is dictated by 
the dust-sound frequency Ud which is the dominant one. It should however be noted that 



UJd is not the usual dust-acoustic frequency which can be written as ^Td/md- The WKB 
coefficients^ a\''^ can now be ordered as as = af\ a2 = 02°"*, cti = ap"* + ai^\ and = aQ'\ 
The second order coefficients are given by 



|e"5„+i(t)| « e"-'Sn{t) , 
e^'Sj.+.it) « 1. 



(46) 
(47) 




(48) 






(49) 




'i^dZdO^d^p- 



(50) 



To the zeroth order [N = 0), the dispersion relation Eq. (1441) now becomes 



(34 + 5T{t)){uoiJoY + 56dZdou;l{t){uoiJo) = 0, 



(51) 



which has two solutions for the condensation and acoustic modes respectively. 



0, 



(52) 




(53) 
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where = Tq ^ and the subscripts 'c' and 's' refer to the condensation and sound modes. 
The first order contribution (A^ = 1) to the condensation mode is given by 



3 

V'lc = --^^p- (54) 



So, the condition for a growing condensation mode is given by the condition Qp < 0, 

+ meT^o—— — < 0, (55) 

^eO OleO V eO / p^o 

which is similar to the condition for a static equihbrium^i^. As can be seen from the above 
condition (l55ll . thermal conductivity has a stabilising effect on temperature perturbation. 
The corresponding critical wave number for the condensation mode is given by^, 

2^ rrieUeO d f Ceo\ 



= -e^T^o-:^^^ ^ , (56) 

beyond which the condensation mode is stabilized. Note that the critical wave number has 
become dependent on the equilibrium electron density, which means that as more and more 
electrons get depleted on the surface of the massive dust particles, perturbation with longer 
wavelengths becomes stable^. This is obvious as the cooling is due to the electrons, the 
depletion of electrons by the dust particles prevents localised regions inside a diffuse structure 
from effectively radiating. The dust can, however radiate at a much longer wavelength to 
which the medium may be opaque^»^. Comprehensive effect of the presence of dust particles 
can only be realised when we consider the effect of dust-charge fluctuation in Section. V. 

The first order contribution to the sound mode can be obtained from the corresponding 
dispersion relation as 

SaPuliipos'tpos + As-ipos) + «2°^ {^o'^osY + a?i^is + = 0, (57) 
which can be solved for ipis, 



3 



2a. 



(0) 
3 



V'ls = — M 4°^(2^e + ^c) - r{u, + uJc) - + (58) 



"2 I "O 



2af 2af)' 



which is real. So the existence of a growing sound mode (overstable) at t = depends on 
whether ipis > at t = 0. 



11 



1 

10-2 
10-^ 

10-" 
10-'° 
10-'" 



1 

- 


1 1 1 


1 1 


—-'TSo.oifr,. 




k = a5fc^_,__— — - 




2k, 








k = 5K - 


1 


1 1 1 


1 1 


1 1 



40 50 



0.0012 
0.001 
0.0008 - 
0.0006 
0.0004 
0.0002 - 


-0.0002 

-0.0004 

-0.0006 

-0.0008 - 

-0.001 





40 50 



Figure 1: Comparison between the WKB and numerically computed solutions of Eq. (j30p for a static 
equilibrium. Relative errors for different wave numbers for the condensation mode are shown in the 
first figure. In the second figure, relative errors for the sound mode for wave numbers /cq = O.lkc 
(dashed) and kQ = (solid) are shown. The parameters Uc = 0.01, (3^^ = 0.5, and Td = uj'^'^ . 

1. Validity of WKB approximation 

At this point, we would like to compare the relative errors between the WKB solutions, 
namely (15^ and ( l58l) of the condensation and the sound modes, and the numerically com- 
puted solutions of Eq. (!30|) . As can be seen from FigUl in the absence of expansion, the 
relative deviation of the WKB solutions from the numerically computed solutions is within 
10~^ except for highly unstable perturbations of very long wavelengths, k <^ kf.. For short 
wavelength perturbation, k > kc, the deviation is down to 10~®. 

In the above analysis, we have assumed that the equilibrium heat-loss function^ Ceo oc 
HeoTll'^ so that the critical wave number kc can be expressed as 



3ac 



1/2 



(59) 



where k^ = kt=o. In the above expression, ac and (3^ denote the ratios of the characteristic 
equilibrium cooling time and the thermal conduction time to the dust-acoustic time at the 
beginning of the expansion (t = 0), 

/3x = (60) 



Or 



UJo 



where (oOco, ^xo^ ~ (^c o equilibrium dust-acoustic frequency at t = 0, 



UJn 



ko \ To' 
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To being the initial equilibrium electron temperature [Teo]t=o- We would like to point out 
that the present form of heat-loss function inherently implies net equilibrium cooling as 
there is no heating source. Any non-zero equilibrium value of £,^0 actually determines the 
net cooling. 

In working out the WKB solutions we have assumed massive dust particles^^ with mass 
of the order of lO^^m//, itih being the proton mass and the dust charge number Zdo = 10^. 
The electron to dust ratio is fixed at 5e = 0.1 i.e. a ninety percent depletion of the electrons 
on the surface of the dust particles. The initial electron temperature is assumed to be 
Teo = 10^ °K at t = 0. The initial equilibrium electron density Ueo ~ lO^cm"^ is taken 
to be of the order of the average electron density in multiple-shell PNe (MSPNe)^'^. The 

electron thermal conductivity at this temperature can be written asM^, 

L84X10-5 _ 

Xo = — ] 7 Tq' ergss K cm , (61) 

log Ac 

where log Ac ~ 30 is the Coulomb logarithm. The domain of instability is defined as 

^^cO < T^o < t^o, (62) 

as thermal conductivity is more effective at short scale length. We have assumed a value 
of ttc = 0.01 and f3^ = 0.5, for which the critical wave number, given by Eq. (l59|) . is ~ 
2.7 X 10^^'^ cm~^ corresponding to a critical wavelength Ac ~ 2.3 x 10^° cm ^ 7.6 x 10^^ pc. 
As the average size of a PNe is of the order of 0.1 pc^^, a perturbation wavelength of the 
order of a few thousand kilometers (~ Ac) is quite small compared to the size of the nebula. 
The characteristic time-scale Td ~ 10'^ years, which is about 10% of the average kinematic 
age of a MSPNe^. 

The equilibrium evolution of electron temperature can be determined from Eq. ffTTl) . which 
in terms of the expansion parameter e can be written as, 

1 - i^^cO - Oi)^ ^ 



Teo(t) = To 



ujcot < 1, (63) 



(l + at)2 

where a = e is the rate of expansion, which is constant in this case. The evolution of Teo(t) 
in the absence of expansion is shown in Figl2l Note that relation fl63l) is valid as long as 
t < oj'^Q. In general, equilibrium cooling should vanish as temperature approaches zero. We 
note that in absence of equilibrium cooling (cjco = 0), the time dependence of the equilibrium 
temperature can be expressed as 

Tco(t)=To£(t)-3(^-i), (64) 
13 



Figure 2: Equilibrium evolution of electron temperature Teo(t) in absence of expansion (first dia- 
gram). The isobaric stability criterion is plotted in the second diagram against k/kc and t. For 
clarity, the negative values of Vtp part of the 0,p — k — t surface are hidden which signifies instability. 
0,p in the figure is normalized by the initial dust-sound frequency wq- 

where 7 is the ratio of specific heat, which reduces to Toe(t)~^ for 7 = 5/3 in our case. 
2. Effect of equilibrium cooling and expansion 

Before we proceed further, a few comments and observations are to be made. If we con- 
sider only a static equilibrium with net cooling, any dependence of the equilibrium quantities 
is due to the evolution of equilibrium temperature. We note that as the stability (or insta- 
bility) condition for the condensation mode i.e. condition (l55ll is valid throughout the whole 
time domain, the expression for critical wave number, given by Eq. (!56l) . is dynamic due to 
the time dependence of the equilibrium quantities. The critical wave number is dynamically 
going to change, effectively determined by the detailed temperature dependence of the heat- 
loss function Ceo- So, if we start off initially with a long wavelength stable perturbation, it 
might become unstable at a later time depending on the instantaneous value of the critical 
wave number kc- As an example, we have shown the evolution of the quantity Qp in time 
with different k/k^ in the second diagram of Fig 121 for the parameter range mentioned later 
in this section. As can be seen from the figure, even if the initial value for fip at t = is 
positive (signifying stability for the condensation mode) for say k ~ 2kc, at a later time, 
around t ^ 25rrf, the quantity falls below zero marking the onset of the instability. Similar 
observations can also be made for the sound mode. For the form of the heat-loss function 
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£eo oc rieoTll'^ , the instability condition fl55|) can be written as, 

+ -3^co<0, (65) 

from which we see that both equihbrium coohng and expansion destabihze the condensa- 
tion mode and an otherwise stable mode ultimately may become unstable due to either 
equilibrium cooling or expansion or both. However, the behavior of the growth rate can 
be ascertained only after a detailed calculation. Apparently, the instability condition for 
the condensation mode does not depend upon the presence of dust particles as it is just a 
balance between the electron conductivity and cooling due to electrons, both of which are 
equally affected by the dust particles. 

As mentioned earlier, the existence of overstable mode depends on the positivity of ipis, 
given by Eq. (l58|) . which can be simplified as, 

V'ls = + 6c^c^ r + ^ Qt^c - T~^^ Se. (66) 

As < 1, we can write the condition of unstable acoustic mode as 

^r-'6e < T (^^Ue + 6cUe) T, (67) 

which, in the limit of vanishing background expansion {u^ = 0) and r ~ 1, becomes T^^^e < 
lOuc and for r <^ 1, reduces to Uc > 2r^^. So, the presence of dust {6e < 1) actually 
destabilizes the acoustic mode and a stable acoustic mode can become unstable due to 
the capturing of the electrons by the massive dust particles. This is due to the decrease 
in thermal conduction by the electrons, as more and more electrons are captured by the 
massive dust particles. The relation, Eq. (l66l) without background expansion (a = 0), can 
be expressed as, 

3^cO fo^ , x\ 3 . _ ,.5,2. 



which shows that equilibrium cooling destabilizes the acoustic mode. 

The immediate effect of equilibrium expansion is that the normal modes of the system 
becomes time dependent, which is analogous to change in frequency of vibration in Doppler 
effects. However the stability of a certain mode may depend on balancing between the 
opposing effects. For example, the instability condition for the acoustic mode in the limit 
of very small equilibrium cooling (cjco ^1), can be expressed as, 

7aT QByk'^Se 



(1 + at) (1 + aty 
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Figure 3: Comparison with the classical case2 for static equilibrium (no expansion) in presence of 
equilibrium temperature evolution for the perturbed dust density for the condensation mode (first 
figure) and wave mode (second figure). Rest of the parameters are same as for Figj2j The dashed 
lines indicate the classical Field's behavior. The values of wave number are shown is at t = 0. 

and we see that background expansion too destabilizes the acoustic mode. Physically, ex- 
pansion nullifies the effect of thermal conductivity as thermal conductivity is most effective 
at short scale lengths, leaving the onset of thermal instability entirely on the isobaric cooling. 

The cumulative effect of equilibrium cooling, background expansion, and electron cap- 
turing my massive dust particles is to destabilize the condensation and the acoustic modes 
which effectively increases (reduces) the cut-off wavenumber (wavelength) beyond which the 
modes are stable. Therefore, we conclude that as the nebulae evolve due to expansion, the 
microstructures eventually die out and the presence of dust particles enhances this effect, 
which actually conforms to the observational studies, which report that the microstructures 
are transient compared to the evolutionary life (due to expansion) of the planetary nebulae^. 

We now compute the WKB solutions of the normalized perturbed dust density n(i{t) for 
the parameter range discussed above. The evolution of the perturbed dust density with only 
equilibrium cooling for the condensation mode and the acoustic mode are shown in Fig|3l 
As can be seen from the figures, though the classical condensation mode is stable beyond 
the critical wave number fcc, in presence of equilibrium cooling, the condensation mode can 
grow further. However, for k ^ fee, the mode is stabilized. The equilibrium cooling has 
similar effects on the acoustic mode (see the second figure in Fig|3]). The classical results 
with no equilibrium expansion and cooling- are shown as dashed curves. The initial level of 
perturbation is fixed at 0.01 in all the cases. 

In FigHl we show the effect of equilibrium expansion on the modes. The rate of equilib- 
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Figure 4: In the first diagram the comparison of the WKB solutions for the condensation mode 
with the classical Field's case in presence of equilibrium expansion and net equilibrium cooling, 
is shown for k = O.Olekc- Other parameters are same as before and the dotted line shows the 
prediction of the classical theory. In the second diagram, the destabilisation of the condensation 
mode for rapid expansion (a = 0.1) is shown for different k. 

rium expansion is given by a. As seen from the first diagram of FigHl rapid expansion causes 
the growth rate of the condensation mode to decrease in the initial phase. However, finally, 
expansion destabilises an otherwise stable mode, which is shown in the second diagram of 
FiglU where we have set the expansion rate at ~ 10% of the dust-acoustic frequency. 



C. The heat-conduction domain (ojc, a;^, „, lj^^ <C t"^^) 

At very small spatial scale, the heat conduction through the electrons dominates and the 
thermal equilibration of the electrons prevents any thermal instability. With this ordering, 
the WKB expansion yields one condensation mode 



Oc 



(3<5e + 5r) 



(70) 



which is stable. 



V. DUST CHARGE FLUCTUATION : NUMERICAL SOLUTION 

In this section we let the dust charge to change according to the dust charging equation 
Eq. (12^ . In what follows, we however neglect the natural decay rate of charge fluctuations 
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Figure 5: Row-wise from the top — (a) Low-frequency response of condensation mode for low 
cooling and expansion, 7c — 7/i = W~'^,a = 10~^ and (b) high cooling and expansion, jc — "fh = 
0.1, a = 0.1. (c) High frequency response to the condensation mode (i) a = 0, 7c — 7h = 0, = 0.1, 
(ii) a = 0.01,7c — 7/i = 0.01, z/ = 0.1, (m) a = 0.1,7c — 7/i = 0.1, = 0.1. (d) Low frequency 
response to the sound mode (i) a = 10~^, 7c — Jh = 10~^, = 0, (ii) a = 10~^, 7c — 7/i = 10~^, u = 
0.01, (m) a = 0.01, 7^ - 7/^ = 0.01, u = 0.01, (iw) a = 0.1, 7^ - 7/^ = 0.1, u = 0.01. 

77. Using Eqs.f ll6p and ( JT5l) we can write the charging equation as 

^ = J^{-rPe-ne), (71) 

where is normalized perturbed dust-charge number and u is the charge fluctuation fre- 
quency normalized to initial dust-acoustic frequency Uq. The fluctuation of dust charge 
introduces another time-scale into the problem, which we would like to keep at an arbi- 
trary level. Therefore a clear separation of time-scales is not possible and we investigate the 
dispersion relation numerically instead of applying WKB approximation. 

In order to facilitate equilibrium balance between radiation loss and heating, we assume 
a form of the net cooling function with both heating and cooling having same dependence 
on density^, 

Ceo = /lKo)[/2(Teo) - ^(Teo)], (72) 
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and the equilibrium temperature Tq is determined by the condition £eo = 0. With the heat- 
loss function given by Eq. fl72|) . the isobaric stability condition becomes /2(Teo) > f^(Teo)- 
Without loss of any generality, we assume that /2 oc T^o, /s oc Tj^, and /i oc Ueo so that the 
net heat-loss function can be written as, 



where Ci,2 are proportionality constants and A and T are equivalent equilibrium cooling and 
heating. The cooling time-scale in terms of Uco can be written as 



with 'fch being the normalised cooling and heating functions and uq = neo\t=o- 

The comprehensive effect of dust-charge fluctuation on the perturbation amplitude is 
shown in FiglS] In the top panel of FigjSl the low-frequency response of the condensation 
mode is shown for low cooling and slow expansion (first diagram) and high cooling and 
rapid expansion (second diagram). Though the effect of charge- fluctuation is to stabilize 
the mode, the effect is more pronounced in the case of rapid expansion as can be seen from 
the level of perturbation amplitude. For charge fluctuation frequency as low as 10% of the 
dust-acoustic frequency, the condensation mode is completely stabilized beyond 20Tfi. As 
background expansion tends to enhance the effect of cooling, we have assumed the cooling 
and expansion time-scales to be of the same order. The high-frequency response of the 
condensation mode is shown in the first diagram of the lower panel of FigjS] where the 
effect of dust-charge fluctuation is shown for different rates of cooling and expansion. We 
observe that with finite dust-charge fluctuation, expansion and equilibrium cooling actually 
destabilise the condensation mode [curve (m) in Figj5]^c)]. However rapid expansion leads 
to stabilisation. This is a distinct departure from the behaviour of the mode in case of 
low-frequency perturbation. 

The behaviour of the acoustic mode with low-frequency perturbation is shown in the 
second figure of the lower panel of FigJSl Increasing equilibrium cooling and expansion 
tend to stabilise the acoustic mode with finite dust-charge fluctuation. High-frequency 
perturbations of the acoustic mode are always damped. 



^eO — n-eoiClTeO — C2T^q) — A — F, 



(73) 




(74) 
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VI. CONCLUSION 



We have considered the the role of thermal instability of an ionized plasma with an 
expanding equilibrium and net equilibrium cooling, which has a considerable presence of dust 
and discussed the outcome of the analysis in the context of expanding planetary nebulae. 
We have found that the cumulative effect of equilibrium cooling, background expansion 
and the electron capturing by the massive dust particles is to destabilize the condensation 
and the acoustic modes which reduces the cut-off wavelengths beyond which the modes 
are stable. Therefore, we see that as a planetary nebula evolves due to expansion, the 
microstructures eventually die out, which is in agreement with the observational reports 
that these microstructures are transient compared to the evolutionary life of a nebula. We 
have shown that the presence of dust particles in this ionized environment tends to destabilize 
the acoustic mode. 

However, the fluctuation of dust-charge severely damps these instabilities and the modes 
are completely stabilized for sufficiently large fluctuation frequency. As there are no conclu- 
sive report on the dust-charge fluctuation parameter in these ionized plasmas, it will be too 
much far fetched to say that the dust-charge fluctuation completely suppresses the formation 
of the microstructures, which can be formed out of radiation induced instabilities. It will be 
of much interest to view these results with newer experimental data. 

Appendix 

In this appendix, we show that for non-static equilibrium, it is essential for the normal 
mode of the system to be expressed exclusively as a function of time and a decomposition^^ 
of the kind expressed in Eq. (12^ is not fruitful. As an example, we take the case represented 
in Fig.l, for k = 0.1 /cc for condensation modes and k = 5kc for acoustic modes. We consider, 
only the effect of equilibrium expansion in the limit of no net equilibrium cooling. Thus the 
equilibrium evolution of the temperature is due to expansion only [see Eq. lITT]) ]. In Fig.6, 
we show the growth of the perturbed dust density nd{t) as obtained by the full numerical 
solution of Eqs.([TlD-([l6D, ([22D, ([24D, ([26]), and ([27D and the solution through the WKB 
expansion, Eq. flHl) . The solution, as given by the prescription represented in Eq. fl2^ is 
plotted in the figures as dashed lines. As can be seen, the solution from the normal mode 
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Figure 6: The numerical solution of Eqs.([T4l)-([T6|), ([22]), ([Ml), dlS]), and dST]), WKB solution, and 
the decomposed solution (shown as dashed lines). The first figure is for condensation mode and 
the second one is for wave mode. The solid curves represent the numerical and WKB solutions 
which are almost indistinguishable. In the same parameter regime, the decomposed solution has a 
large deviation in amplitudes and phases. 



decomposition has a large deviation from the exact numerical solution, in terms of growth 
rates, whereas, the WKB solution almost coincides with the numerical solution. Also note 
the difference in phases of the wave mode solution from the numerical one, while the WKB 
solution has been able to keep in phase with the numerical solution. 
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